MHD Energy Bypass Scramjet 
Performance With Real Gas Effects 


Chul Park* 

ELORET Corporation, Sunnyvale, CA 94087 
Unmeel B. Mehta* 

NASA Ames Research Center, Moffett Field, CA 94035 

and 

David W. Bogdanoff x 

ELORET Corporation, Sunnyvale, CA 94087 ‘ 


Abstract 


Copyright @ 2000 by American Institute of Aeronautics and Astronautics, Inc. No copyright is asserted in the 
rio? a un r?H Under ! 1 * U S - COde - The U S - Government has a royalty-free" license to exercise a^ 
copyrigW^ner 6 C ° Pyr ' 9h C a,med herein ,or government Purposes. All other rights are reseived by the 

Senior Research Scientist; Fellow, AIAA; e-mail: cpark@mail.arc.nasa.gov; mailing address: Mail Stop 
229-1 , NASA Ames Research Center, Moffett Field, CA 94035. 

Division Scientist, Space Technology Division; Associate Fellow, AIAA; e-mail: umehta@mail.arc.nasa.gov; 
mailing address: Mail Stop T27B, NASA Ames Research Center, Moffett Field, CA 94035. 

Senior Research Scientist; Associate Fellow, AIAA; e-mail: dbogdanoff@mail.arc.nasa.gov; mailing 
address: Mail Stop 230-2, NASA Ames Research Center, Moffett Field. CA 94035. 


The theoretical performance of a scramjet propulsion system incorporating an magneto-hydro-dynamic 
(MHD) energy bypass scheme is calculated. The one-dimensional analysis developed earlier, in which the 
theoretical performance is calculated neglecting skin friction and using a sudden-freezing approximation for 
the nozzle flow, is modified to incorporate the method of Van Driest for turbulent skin friction and a finite-rate 
chemistry calculation in the nozzle. Unlike in the earlier design, in which four ramp compressions occurred in 
the pitch plane, in the present design the first two ramp compressions occur in the pitch plane and the next 
two compressions occur in the yaw plane. The results for the simplified design of a spaceliner show that (1) 
the present design produces higher specific impulses than the earlier design. (2) skin friction substantially 
reduces thrust and specific impulse, and (3) the specific impulse of the MHD-bypass system is still better 
than the non-MHD system and typical rocket over a narrow region of flight speeds and design parameters. 
Results suggest that the energy management with MHD principles offers the possibility of improving the 
performance of the scramjet. The technical issues needing further studies are identified. 


Nomenclature 


B: Magnetic field strength, Tesla. 

E: Electric field strength, V/m. 

H: Height of vehicle, m. 

h: Reactor entrance height/width, m. 

M e : Combustor entrance Mach number 

P: Pitch (Fig. 3). 

T: Temperature, K. 

x: Axial distance, m. 

y: Lateral distance, m. 

u: Axial velocity, m/s. 

V: Flight velocity, m/s. 

a,: Load factor (E/uB) for generator. 
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a 2 : Load factor (E/uB) for accelerator. 

5: Boundary layer momentum thickness, m. 

9: Ramp angle, deg. 


Introduction 


In a scramjet propulsion system, one well-known difficulty is the fact that combustion is inefficient in the 
combustor when the Mach number of the airflow entering the combustor is high. In Refs. 1-3. it was noted 
that energy management with magneto-hydro-dynamic (MHD) techniques present a possible means for 
extending the flight Mach number envelope of conventional engines. This idea can also be used to improve 
the performance of the scramjet combustor by the use of the MHD principles. 4 In the proposal, the flow 
entering the combustor is to be decelerated by the use of an MHD generator. The electrical energy extracted 

from the MHD generator is expended to accelerate the flow by the use of an MHD accelerator after the * 
combustion and prior to entering the nozzle expansion. 

In Refs. 4 and 5. the theoretical performance of such an MHD-energy bypass scramjet propulsion system 
is calculated under several simplifying assumptions. The most important of those assumptions are that: 1) 
the flow is inviscid and therefore no friction drag is produced. 2) the chemical reactions in the expanding 
nozzle undergo freezing suddenly, and 3) inlet compression occurs by four ramp compressions in the pitch 
plane. The calculation presents an envelope of specific impulse values attainable with such a system under 
these assumptions, and compares it with those for a non-MHD scramjet system. The comparison shows 
that, in the flight speed range from about 3.4 to 4 km/s. the MHD-bypass system produces specific impulses 
higher than the non-MHD system and a typical rocket. 

In the present work, the earlier method of performance calculation is improved. The flow is taken to be 
viscous, and the chemical reactions occur at a finite rate. The first two ramp compressions occur in the pitch 
plane, but the next two ramp compressions occur in the yaw plane. Other assumptions used in Ref. 4, such 
as ideal one-dimensional performance of the MHD devices, are retained in the present work. 
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The results show that finite-rate chemistry produces nearly the same results as the sudden-freezing 
approximation, the new geometry produces generally larger specific impulse, but the friction substantially 
reduces specific impulse. However, compared with the non-MHD system, the MHD-bypass system is still 
superior at a speed range above 3.4 km/s under certain restrictions on the geometry. The areas needing 
further research and development are identified. 

Method 

Two-Plane Four-ramp Design 

In an idealized scramjet system, the compression of the air flow entering the combustion chamber is 
achieved through the use of a curved ramp, which compresses air through a nearly isentropic process to an 
. an9 ' e t yP' cal, y of about 15 degrees. At the end of this compression, the flow is turned abruptly in the inlet of 
the engine so that the flow becomes approximately parallel to the oncoming free stream. The flow then 
enters an isolator wherein the flow passes through multiple weak shock waves. The isolator isolates the 
effect of combustion-induced disturbances on inlet operation. The flow is generally three-dimensional. The 
combustor following the isolator is usually relatively long, because a long distance is needed in order for the 
flow to complete combustion at high supersonic or hypersonic speeds. 

In Ref. 4, this complex compression process is approximately represented by a simple four-ramp system. 

In that representation, the flow makes four turns of a same specified angle within the pitch plane, as shown 
schematically in Fig. 1. This design will be referred to hereafter as single-plane four-ramp compression 
design. In this design, the flow entering the reactor (MHD generator + combustor + MHD accelerator) has a 
laterally elongated cross section: the width spans the entire vehicle and the height is only about few 
centimeters. In an inviscid calculation, this elongated cross-section had no consequences. 
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Figure 1. Quasi-one-dimensional single-plane 
four-ramp compression inlet system. 4 

U» skin friction inside such an elongated reactor becomes large because the geometry produces a large 
wetted surfaced area. In order to reduce the skin friction to the minimum, the wetted area Inside the reactor 
must be kept to the minimum. The minimum wetted area is obtainable with a square cross section. When 
other design factors are considered, a square cross section could be non -optimum or inappropriate. 



Figure 2. A schematic of two-plane four-ramp 
compression Inlet system with shock patterns. 

As In the single-plane lour-ramp design, all turn angles are assumed to be equal. The consequences of 
deviating from this assumption are discussed later. The first two turns are made within the pitch plane, but 
the next two turns are made within the yaw plane, as shown schematically in Figs. 2 and 3. The lengths ol 
ramps are specified in Ms model such that, at the end ol the lour-step compression, the cross-section ol the 
flow path is a square. The dimensions ol the design are uniquely determined by specifying the common 
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ramp angle 9 and the height of the vehicle H (see Fig. 3). The overall length of the vehicle is fixed at 46 m. 
Of particular interest is the pitch, P, and the width and height of the reactor at its entrance, h (see Fig. 3). 
This design is termed two-plane four-ramp compression design. The scram, et vehicle employing such a 
four-ramp system is shown schematically in Fig. 4. Several of those engine units are ganged laterally. 

In Fig. 5. the reactor is shown schematically for both non-MHD and MHD-bypass systems. For a non- 
MHD system, the reactor consists of an isolator, a fuel injector, and a combustor. For an MHD-energy 
bypass system, the isolator must function also as a seeder and an MHD power generator. A fuel injector 
must follow the MHD generator. The combustor, connected to the fuel injector, must function also as an 
MHD accelerator. The lengths of these components are chosen arbitrarily as shown in Fig. 5 for the present 
work. The cross-section is assumed to be always a square. The dimensions of the square vary slightly for 

the MHD case as shown schematically, because of the increase in entropy due to Joule heating and 
because of the flow acceleration. 




Figure 3. Flow path in the two-plane 
four-ramp compression Inlet system. 
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Figure 4. Schematic view of a scramjet vehicle 
employing two-plane four-ramp inlet system. 


Inviscid Calculation Along Flow Path 


As in Ref. 4, the flow properties are assumed to be in equilibrium up until the point in the nozzle at which 
chemical freezing takes place. In the MHD generator and accelerator, the flow is assumed to be uniform 
across the cross section; and the ideal one-dimensional MHD equations, given in Ref. 4, are assumed to be 
valid. The parameter dictating the design of the MHD devices is the load factor a = E/uB. This parameter, 
the dimension of the entrance of the MHD device (h for the generator and the corresponding dimension for 
the accelerator, see Fig. 3). and the requirement that the MHD action be completed at the end of the given 
length of the MHD device completely specify the design of the MHD devices including the magnetic field 

strength B. electrical field strength in the lateral and axial directions. E, and E,. electrical current, and the 
Hall parameter. 

At the end of the fourth compression, seeding is assumed to occur. Liquid potassium or cesium is 
assumed to be injected into the flow instantly without causing shock waves or viscous dissipation. The flow 
exiting from the MHD generator enters the combustion chamber. Here, gaseous hydrogen fuel is assumed 
to be mixed instantly without causing any shock waves or viscous dissipation. 

In Ref. 4, the mass fraction of the seed substance was varied. The specific impulse was found to be only 
mildly dependent on the seed mass fraction. The optimum values of mass fraction were found to be about 

0.001 for potassium and 0.003 for cesium by mass. These values are assumed in the present work 
throughout. The method of fuel injection is described in Ref. 4. 
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For the purpose of calculating the inviscid flow properties, the expansion following the MHD accelerator 
is assumed to occur with a linearly increasing cross-sectional area. The rate of increase of the cross- 
sectional area is chosen so that, at the end of the expansion, the height of the nozzle equals the height of 
the vehicle, as indicated in Fig. 4. 
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Figure 5. Schematic of the reactor. 

In Ref. 4, the flow properties are calculated for the expansion region using a sudden-freezing 
approximation. In the present work, this is replaced by a finite-rate kinetic calculation. The chemical species 
considered are O, N, H, 0 2 , N 2 , H 2 , NO, OH, H 2 0, and a seed substance, its ions, and electrons. The seed 
substance is potassium, K, or cesium, Cs. The forward (endothermic) reaction rates are expressed as 

k = AT e' TrfT cm 3 /(mol-s). 

The parameters A, n, and Tr are obtained from Refs. 6 to 9. They are summarized in Table 1. 

Table 1. Reaction rates used. A is in the units of cm s /(mol-s). M is the third body. 
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M 

A n 

Tr 

Reference 

) H 2 + M -» H + H+M 


0 

2.940e 1 8 -1.0 

52500 

Baulch ©t al 7 

N 

2.940e 18 -1.0 

52500 

Baulch et al 7 

H 

2.940e 18 -1.0 

52500 

Baulch et al 7 

K/Cs 

2.940e 18 -1.0 

52500 

Baulch et al 7 

0 2 

2.940e 18 -1.0 

52500 

Baulch et al 7 

N2 

2.940© 18 -1.0 

52500 

Baulch et al 7 

h 2 

7.350© 18 -1.0 

52500 

Baulch et al 7 

NO 

2.940© 18 -1.0 

52500 

Baulch et al 7 

OH 

2.940© 18 -1.0 

52500 

Baulch ©t al 7 

H 2 0 

4.780© 19 -1.0 

52500 

NASP cmt.® 

K7Cs 

** 2.940e 18 -1.0 

52500 

Baulch et al 7 

1 0 2 + 

M —> 0 + 0 + 

M 


0 

5.090© 18 -1.1 

59360 

Park® 

N 

5.090© 18 -1.1 

59360 

Park 6 

H 

5.090© 18 -1.1 

59360 

Park 6 

K/Cs 

5.090© 18 -1.1 

59360 

Park 6 

0 2 

5.090© 18 -1.1 

59360 

Park 6 

N2 

5.090© 18 -1.1 

59360 

Park 6 

h 2 

2.000© 18 -1.1 

59360 

Park 6 

NO 

1.270© 19 -1.1 

59360 

Park 6 

OH 

5.090© 18 -1.1 

59360 

Park 6 

H 2 0 

8.270© 19 -1.1 

59360 

Park 6 
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KVCs 

* 5.090e 18 -1.1 59360 Park 5 

(3)OH + M -» + H + M 


0 

2.400© 13 0.3 5000 

Baulch et al 7 

N 

2.400e 13 0.3 5000 

Baulch et al 7 

H 

2.400e 13 0.3 5000 

Baulch et al 7 

K/Cs 

2.400e 13 0.3 5000 

Baulch et al 7 

0 2 

2.400e 13 0.3 5000 

Baulch et al 7 

N2 

2.400e 1 3 0.3 5000 

Baulch et al 7 

h 2 

6.000e 13 0.3 5000 

Baulch et al 7 

NO 

2.400e 13 0.3 5000 

Baulch et al 7 

OH 

2.400e 1 3 0.3 5000 

Baulch et al 7 

HjO 

3.900e 14 0.3 5000 

Baulch et al 7 

KVCs* 

2.400e 13 0.3 5000 

Baulch et al 7 

(4) H 2 0 + 

M OH + H + M 


O 

1 .060e 25 -2.5 61000 

Baulch et al 7 

N 

1 .060e 25 -2.5 61000 

Baulch et ai 7 

H 

1 .060e 25 -2.5 61000 

Baulch et al 7 

K/Cs 

1 .060e 25 -2.5 61000 

Baulch et al 7 

0 2 

1.060e 25 -2.5 61000 

Baulch et al 7 

N2 

1 .060e 25 -2.5 61000 

Baulch et al 7 

h 2 

2.650e 25 -2.5 61000 

Baulch et al 7 

NO 

1 .060e 25 -2.5 61000 

Baulch et al 7 

OH 

1 .060e 25 -2.5 61000 

Baulch et al 7 

O 

N 

X 

1.720e 26 -2.5 61000 

Baulch et al 7 

KVCs* 

1 .060e 25 -2.5 61000 

Baulch et al 7 

(4) K 

+ e~ -* K* + e' + e' 
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3.900e 33 -3.78 45180 Park 6 

(4) Os + © — > Cs* + 0 + 0 ~ 

3.900©33 -3.78 45180 Park 6 

(5) N 2 + 0 -> NO + N 

5.6900 12 0.42 42938 Bos© 0 t al 9 

(6) O + NO -» N + 0 2 

2.360© 9 1.00 19220 Park 6 

(7) H 2 0+ H H 2 + OH 

2.760e 10 1.12 10250 Baulch e t al 7 

(8) H 2 + 0 H + OH 

5.060© 4 2.67 3166 Baulch etal 7 

(8) 0 2 + H -> OH + O 

LS10e 14 0, 8273 Baulch fit al 7 

Skin Friction 

Tha boundary layer flow along th© wall in th© flow path is assumed to be fully turbulsnt. Van Driest 
equations’ 0 " for fully turbulent flow over a flat plat© are used to calculate the skin friction. These equations 
yield the skin friction over a flat plate between the given distances x, and x 2 away from the sharp leading 
edge. At x„ the equations yield also the boundary layer momentum thickness 5,. The thickness 8, can be 
used interchangeably to denote the distance x t . This relationship between 8,. rather than x„ and skin friction 
is made into a routine in the present work. 

This routine is then used in calculating the friction over a curved surface or a flat surface with varying flow 
conditions. The surface is divided first into small segments of length Ax = x 2 - x, = 0.05 m. At the apex of 
the vehicle or the leading edge of the cowl (i.e. the second ramp, see Fig. 3). 8 is zero. Over each Ax, the 
flow properties at the edge of the boundary layer are assumed to be unchanged. Then, using the routine, the 
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friction over that segment is calculated. The calculation yields the 5 value at the exit of this segment, S 2 . This 
5 2 value is used in the calculation for the next segment, and so on. 

In the region of the third ramp (the hatched area in the bottom figure in Fig. 3), the streamline starting at 
different lateral positions encounter the third and the fourth shock at different locations. For this region, the 
friction over a unit width is first calculated along five different streamlines, and is integrated to evaluate the 
total friction over these panels using the Simpson’s rule. Calculations showed that five streamlines lead to 
solutions that are accurate to within 0.1% compared to those solutions with nine streamlines. 

For the contoured portion of the expansion, shown schematically in Fig. 6, accurate calculation of skin 
friction is difficult because the contours are unknown. Here, skin friction is assumed to be that of the four flat 
plates shown in Fig. 7. The lengths of these flat plates are arbitrarily assumed to be 3 m as shown. The 
inaccuracies in computing the skin friction drag do not impact the conclusion regarding which design, MHD 
or non-MHD, generates higher specific impulse, since drag in both cases is computed in the same manner. 



Figure 6. Schematic of the transition from 
the MHD accelerator to one-sided expansion. 
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Figure 7. Flat plates assumed for friction calculation of the contoured expansion region. 


Highest Specific Impulse 

The combustor entrance Mach number M e is varied as 1.5, 2.0. and 2.5, which are within the realistic 
range of possibilities for a scramjet combustor. The vehicle height H is varied as 3.5, 4, 4.5, 5. and 5.5 m. 
although most of the results are for 3.5 m. The load factor pair, a, - a 2 - is varied as 0.95-1/0.95, 0.90-1/0.90, 

and 0.80-1/0.80. The fuel equivalence ratio is kept at unity for all cases, and the flight dynamic pressure is 
kept at 1 atm. 

For a flight speed V, the combustor entrance Mach number M Cf the vehicle height H, and the seed 
substance and seed mass fraction chosen, one can vary the ramp angle 9 arbitrarily. The resulting specific 
impulse value is a function of the ramp angle. Calculation is performed at an interval in 0of 0.1° to first 
produce a table of specific impulse versus 9. The ramp angle giving the highest specific impulse is chosen 
from this table as the optimum ramp angle for that condition. 

The thrust is calculated as the difference between the outgoing (at the tail of the vehicle) and incoming 
momenta (at the nose of the vehicle). Specific impulse is calculated by dividing this thrust by the sum of the 
flow rates of the fuel and the seed material. 
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Results 


General Features of Solutions 


in Fig. 8. the specific impulses obtained by the present method are compared with those obtained by the 

sudden-freezing approximation used in Refs. 4 and 5. for the singfe-pfane four-, amp compression system 

considered therein. The combustor entrance Mach number M c was 2 for this case. Skin friction drag Is no. 

incuded in mis comparison. As seen in Fig. 8. there is only a small difference between the kinetic and 
sudden-freezing calculations. 



Figure a. Comparison of the specific Impulse obtained by kinetic calculation and that with the 
sudden-freezing calculation; or, = 0.95, a, = 1/0.95; and potassium seed. 
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Compared with the single-plane four-ramp design of Ref. 4. the present two-plane four-ramp design 
generates a vehicle of a smaller height, H. As a result, the area ratio of nozzle expansion is slightly smaller 
than in the previous design. The length of the nozzle is affected by the height: the larger the height, the 
shorter becomes the nozzle. For the vehicle of H less than 5 m, the present nozzle is longer than that of the 
single-plane four-ramp design. As a result, for the present vehicles of H < 5 m, expansion becomes more 
gradual than for the single-plane four-ramp design. This gradual expansion results in a nozzle flow closer in 
chemical equilibrium at the exit. This in turn produces a greater thrust. However, the longer nozzle also 

produces greater friction drag. Thus the functional relationship between the height H and specific impulse 
becomes somewhat complicated. 

In Fig. 9, the optimum vehicle geometry is compared between typical MHD and non-MHD cases with 
friction and with finite-rate chemistry in the nozzle. For this case, V = 3.75 km/s, a, and a 2 are 0.95 and 

1/0.95, M c = 1.5, and H = 3.5 m. As indicated, the optimum angle of the ramps for the MHD and the non- 
MHD cases are 9= 16.5 and 25.6°, respectively. As has been pointed out in Ref. 4, this difference in the 
optimum 0 occurs because of the constraint on M c . 

The pressure distribution along the flow path for these two vehicles is compared in Fig. 10. As shown 
here, because of the large 0, the non-MHD vehicle is subject to a higher pressure in the inlet region than the 
MHD vehicle. This causes larger form drag for the non-MHD vehicle. The specific impulses for these two 
cases are 791.7 sec for the MHD vehicle and 704.7 sec for the non-MHD vehicle. 





Figure 9. Comparison of the optimum configuration between the MHD and the non-MHD vehicles. 
V = 3.75 km/s, M e = 1.5, H = 3.5 m, a, = 0.95, a, = 1/0.95, potassium seed. 


The friction drag per meter of the vehicle width to the two vehicles are shown by components in Figs. 
11(a) and (b). Because 0 is different between the two vehicles, the friction drag values are slightly different 
on the three ramps. The friction drag for these ramps is slightly larger for the MHD case because the MHD 
case produces longer ramps. The sum of the drags by the MHD generator, fuel injector, and the MHD 

accelerator for the MHD vehicle, shown in Fig. 1 1 (b), is lower than the sum of the drags by the isolator and 
the combustor for the non-MHD vehicle shown in Fig. 1 1 (a). 

The details of a typical solution for a two-plane four-ramp compression design are given in Table 2. The 
highlights ol mis sol ut ion are .ha.: (,) ,he heigh, is 3.5 m. «he ramp angte is ,6.5-, and ,he pitch is ,.842 m; 
(2) me vehic,, produces 43 tons ol thrus, per , meter o, vehicle widm and a specilic impulse o, 791.7 sec 
(as opposed to 704.7 sec lor the non-MHD vehicle); (3) the load lactont a, and ct,are 0.95 and 1/0.95, 
respectively. (4) me magnetic fields required lor the MHD generator and accelerator are 11.29 and 9.115 

Tesla, respectively; and (5) the axial voltage gradient, which must be below about 5000 V/m (see Hel. 4). is 
below 4425 V/m. 
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Pressure, Pascal 



Flour. 10. Comparison ol static pressor, distribution ov.r the MHO and the non-MHD vehicles with 
friction and finite-rat. chemistry m th. noalo. V =. 3.75 km/s, M,= 1.5, H . 3.5 m, o, = 0.95, a,. 1/0.95 




Figure 1 1 . Breakdown of friction drag per 1 m of vehicle width; 
V = 3750 m/s,M e = 1 .5, H = 3.5 m. (a) non-MHD vehicle. 



0 5000 1 10 4 1.5 10 4 2 10 4 2.5 10 4 3 10 4 3.5 10 4 4 10 4 

Friction drag par mater of width, Newton 


Figure 11(b). MHD vehicle; a, = 0.95, a,= 1/0.95; potassium seed. 
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Table 2. Summary of a typical solution. 


Freestream 
Flight speed 
Flight dynamic pressure 

Mach number 
Freestream temperature 
Freestream density 
Freestream pressure 
Number of ramps 
Ramp angles 
Over 1st ramp 
Pressure 
Temperature 
Velocity 
Mach number 
Pitch 

Over 2nd ramp 
Pressure 
Temperature 
Velocity 
Mach number 
Over 3rd ramp 
Pressure 


3.75 km/s 
1 atm (2116 psf) 
(1.014x10* Pascal) 
11.81 
250 K 

1.441 x10' z kg/m 3 
1. 039 xIO 3 Pascal 
4 

1 6*5 * 


2.303 xIO 4 Pascal 
1378 K 
3487 m/s 
4.674 
1.842 m 

1 .095 x 1 0* Pascal 
2145 K 
3190 m/s 
4.048 


3.954x10* Pascal 
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Temperature 


2848 K 


Velocity 3154 m/s 

Mach number 3.451 

i 

Over 4th ramp (=entrance to seeder) 

Pressure 1.113x10® Pascal i 

Temperature 3397 K 

Velocity 2792 m/s 

Mach number 2.768 

Capture air flow rate 1 89. 1 5 kg/s 

Enthalpy flow 1.565x10® W/m 

Seeding 

Entrance height/width 33.50 cm 

Overall inlet area ratio 51 .96 

Seeding material potassium 

f 

Seed mass fraction 10' 3 I 

Seed flow rate 0.189 kg/s 

Pressure 1.116x10® Pascal 

Temperature 3397 K 

Ionization mol fraction 4.431x10'* 

Electrical conductivity 33.25 mho/m 

Entrance to MHD generator 
Height/Width of channel 33.50 cm 

Magnetic field 1 1 .29 T esla 

Hall parameter 2.808 

T ransverse voltage gradient -29,920 V/m 
Axial voltage gradient 4421 V/m 

Voltage across electodes 1 0020 V , 
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Current density 

-5.236 xIO 4 A/m 2 

Pressure 

1.116x10® Pascal 

Temperature 

3397 K 

Velocity 

2789 m/s 

Mach number 

2.764 

Enthalpy flow 

1.676 x10® W/m 

At exit of MHD generator 

Magnetic field 

11.29 Tesla 

Hall parameter 

2.970 

Transverse voltage gradient 

-16,140 V/m 

Axial voltage gradient 

2523 V/m 

Voltage across electrodes 

7572 V 

Current density 

-2.876x10® A/m 2 

Height/Width o? channel 

46.91 cm 

Pressure 

1.054x10® Pascal 

Temperature 

3393 K 

Velocity 

1 505 m/s 

Mach number 

1 492 

Length of generator 

2.721 m 

Enthalpy flow 

1.043x10® W/m 

Power extracted 

5.220x10® W/m 

At entrance to combustor 

Equivalence ratio 

1 

Fuel flow rate 

5.553 kg/s 

Height/Width of channel 

46.91 cm 

Pressure 

1.054x10® Pascal 

Temperature 

3393 K 
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Velocity 
Mach number 

Fuel injector 
Fuel 

Total temperature 
Total pressure 
Nozzle exit velocity 
Nozzle exit static pressure 

At exit of injector section 
Pressure 
Temperature 
Velocity 
Mach number 
Enthalpy flow 

At entrance to MHD accelerator 
Magnetic field 
Hall parameter 
Transverse voltage gradient 
Axial voltage gradient 
Voltage across electrodes 
Current density 
Height/Width of channel 
Pressure 
Temperature 
Velocity 
Mach number 
Ionization fraction 


1 505 m/s 
1.492 


gaseous hydrogen 
500 K 

2x10 7 Pascal 
2865 m/s 
1.054x10® Pascal 


1.035x10® Pascal 
3518 K 
1 555 m/s 
1.327 

1.082x10® W/m 


9.1 15 Tesla 
2.212 

14,920 V/m 
1650 V/m 
8129 V 

2.597 xIO 4 A/m 2 
54.49 cm 

1.035x10® Pascal 
351 8 K 
1 555 m/s 
1.327 

5.302x10* 
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Electrical conductivity 

34.82 mho/m 

Enthalpy flow 

1.080 x10 s W/m 

At exit of MHD accelerator (=entrance to nozzle) 

Magnetic field 

9.1 15 Tesla 

Hall parameter 

1.317 

Transverse voltage gradient 

26,970 V/m 

Axial voltage gradient 

3124 V/m 

Voltage across electrodes 

11,190 V 

Current density 

4.724 xIO 4 A/m 2 

Height/Width of channel 

41 .49 cm 

Pressure 

9.859 x10 s Pascal 

Temperature 

3502 K 

Velocity 

281 1 m/s 

Mach number 

2.400 

Length of accelerator 

2.846 m 

Enthalpy flow 

1.612x10 s W/m 

Power consumed 

5.334x10® W/m 

At nozzle exit 

Area ratio 

37.7 

Pressure 

9. 173x10® Pascal 

Temperature 

1951 K 

Velocity 

4124 m/s 

Enthalpy flow 

1.604 x10 s W/m 

Mach number 

4.441 

Overall performance 

Energy bypass ratio 

0.3336 

Net thrust 

4.31 IxlO 4 N 


23 


Seed + fuel flow rate 


5.553 kg/s 
Specific impulse 79 1 .7 sec 


Parametric Study 

In Fig. 12, the calculated specific impulse values for M c = 1.5 are compared between the MHD case with 
potassium seeding and cesium seeding and with a, =0.95 and a 2 = 1/0.95, the non-MHD case, and a typical 
rocket (with specific impulse of 450 sec). The plotted “with friction" values are obtained by choosing the 
ramp-angle 0 that maximizes the specific impulse while including friction in the thrust calculation. That is, the 
design is optimized including friction. The “without friction" values plotted are obtained by removing friction 
from the thrust calculation for the same vehicle at the same flight condition. The “without friction" values are 
therefore not of the optimized vehicles. 

As seen in Fig. 12, the skin friction for fully turbulent flow greatly reduces the specific impulse. The two 
MHD cases produce higher specific impulses than the non-MHD case between the flight speeds of 3400 and 
4000 m/s. Potassium seeding gives about the same specific impulse as cesium seeding. Compared with the 
single-plane four-ramp results shown in Fig. 8, the present results show higher specific impulses when 
friction is ignored (“without friction" values). This is a result of the gradual expansion in the nozzle mentioned 
earlier. 
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3000 



Figure 12. Comparison of specific impulse with friction between the MHD case with potassium 


and cesium seed and the non-MHD case. M e = 1.5, a, = 0.95, a*= 1/0.95, and H = 3.5 m. 


The effect of changing the combustor entrance Mach number M e on specific impulse is studied in Fig. 1 3. 
As seen here, the case of M c » 1 .5 with MHD yields higher specific impulses than the non-MHD case at flight 
speeds above 3400 m/s. At M e greater than 1 .5, the MHD cases give specific impulses equal to or lower 
than the non-MHD case. 

The effect of the load factors a, and a 2 on specific impulse is studied in Figs. 14(a) to (c). Specific impulse 
deteriorates with a, and a 2 deviating farther away from unity. 

The impact of varying vehicle height H on specific impulse is shown in Fig. 15 for V = 3.75 km/s and M e = 
1.5. By increasing H, the specific impulse of the MHD-bypass vehicle is increased. The specific impulse for 
the non-MHD vehicle also increases slightly more. As a result, for this condition the advantage of the MHD- 
bypass system decreases at large H. 
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ai 
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typical rocket 



3400 
Flight velocity, m/s 


Figure 14. Effect of load factors a, and a 2 on specific impulse. \a) M e = 1.5. 
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Figure 15. Impact of vehicle height H on specific impulse. 


The optimum geometry of the present two-plane four-ramp compression design is one in which the 
reactor is located close to the leading edge, as shown In Fig. 9. The location of the reactor can be moved aft 
by increasing the ramp angle of the third and fourth ramps. Calculation was performed with the angle of the 
third (and fourth) ramp that is 1 .5, 2, and 2.5 times the angle of the first (and the second) ramp. The overall 
length of the vehicle is kept the same. In Fig. 16, the result of that calculation is shown. As the figure shows, 
increasing the third/fourth ramp angle can increase the distance from the vehicle leading edge to the reactor 
entrance. However, specific impulse is decreased significantly due to increased total pressure losses at the 
two stronger shocks. 


Discussion 
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In Ref. 4, a method was presented for assessing the theoretical performance of an MHD-energy bypass 
scramjet propulsion system under several simplifying assumptions. In the present work, this method of 
performance calculation is improved by eliminating two key assumptions. The flow is taken to be viscous, 
and the chemical reactions occur at a finite rate. The improved method provides a more realistic 
assessment of this propulsion concept than the previous method did. 

In the present work, the boundary layer in the flow path is assumed to be fully turbulent starting from the 
leading edge. The boundary layer over the first ramp and the second ramp are no doubt at least partly 
laminar. Consideration of laminar boundary layer over these two components will lower the skin friction 
values for these two components. However, according to Figs. 11(a) and (b). the contribution of these two 
components to the total drag is relatively small. The end result is that the inclusion of viscous friction lowers 
specific impulse substantially. 

With the present design, the MHD-bypass system produces specific impulses that are better than those 
for the non-MHD scramjet system and for a typical rocket engine between the flight velocity of 3400 and 
4000 m/s. Over a narrow range of design parameters, namely, for combustor entrance Mach number of 1.5, 
the vehicle height smaller than 4 m for a 46 m long vehicle, the MHD-bypass system has an advantage over 
the non-MHD system. At the flight velocity of 3900 m/s with M c = 1.5, the specific impulse for the MHD case 
exceeds that of the non-MHD case ant that of the typical rocket by about 150 secs. At this condition, the 
MHD-bypass system offers a clear advantage over the non-MHD case for hypersonic cruiser application. 
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Figure 16. Effect of changing the ratio of the third and fourth ramp angles to 
the first and second ramp angles. 


As seen above, the present two-plane four-ramp compression design is superior to the single-plane four- 
ramp design in Ref. 4. With the present design, performance is seen to be a complex function of the ramp 
angles and vehicle height. In reality, portion of compression will be by an isentropic compression instead of 
shock compression. When Isentropic compression is introduced into the design, performance will be an 
even more complex function of geometry. This offers hopes for further improving the performance of the 
MHD-energy bypass scramjet system. 

In reality, it is probably rather difficult to realize some of the assumptions made for the MHD-system. 
First, the computed results for magnetic field in the magnetic devices are quite high. These values are 
currently not feasible for the spaceliner application. Second, the skin friction through the MHD devices may 
be larger than that for the non-MHD system, because the MHD devices have electrodes which produces 
effectively a rough surface. Third, equilibrium ionization with seeding requires temperatures greater than 
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3000° K, presenting a challenge for designing an appropriate thermal management system and for 
developing appropriate thermal protection system including ultra-high-temperature and light-weight 
materials. Fourth, three-dimensional simulations are essential for fully accounting the effects of the electric 
and magnetic fields on the propulsive flow field. Fifth, specific impulse of airbreathing hypersonic propulsion 
system is very sensitive to vehicle design and operation. 

At moderate hypersonic flight Mach numbers, the energy management with MHD principles facilitates 
operation of scramjet as if it is operated at low hypersonic flight Mach numbers. A great deal of research and 
development would seem necessary for the MHD devices to function as efficiently as assumed. The viability 
of energy management with MHD devices needs to be determined by conducting conceptual designs of 
realistic spaceliners with MHD energy bypass scramjets and assessments of vehicle performance using 
analysis methods with further improvements. 


Conclusions 

The energy management with MHD offers the possibility of enhancing the performance of scramjet. An 
MHD-energy bypass scramjet with a two-plane four-ramp inlet geometry has a higher thrust performance 
than one with a single-plane four-ramp inlet geometry. If the boundary layer of the flow in the propulsive 
flowpath for the MHD-energy bypass scramjet propulsion system is fully turbulent, if the length of the reactor 
for the MHD-system (MHD generator + fuel injector + MHD accelerator) is the same as that of the non-MHD 
scramjet system, and when the height-to-length ratio of the vehicle is limited to below a certain value, then 
the specific impulse of the MHD system is slightly greater than that of the non-MHD system and a typical 
rocket within a speed range between 3400 and 4000 m/s. The specific impulse of the MHD system is nearly 
150 sec higher than that of the non-MHD system and that of a typical rocket systems at 3900 m/s. At this 
speed, the MHD system offers an advantage over the non-MHD system for the hypersonic cruiser 
application. At flight speeds higher than 4000 m/s, the MHD-system is better than the non-MHD system, but 
the absolute value of specific impulse is below that of a typical rocket. 
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Specific impulse is a complex function of geometry; by varying geometry and by improving the analysis 
method, further improvements in performance may be possible. A system analysis of a realistic spaceliner 
with MHD-bypass scramjet is recommended to assess the viability of this propulsion concept. Major 
research and technology issues are identified. 
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